Effect of the shape of periodic forces and second periodic forces on 
horseshoe chaos in Duffing oscillator 

V. Ravichandran a , V. Chinnathambi a , S. Rajasekar b '*, Choy Heng Lai c 

^Department of Physics, Sri K.G.S. Arts College, Srivaikuntam 628 619, Tamilnadu, India 
b School of Physics, Bharathidasan University, Tiruchirappalli 620 024, Tamil Nadu, India 
c Department of Physics, National University of Singapore, Singapore 117542 



Abstract 

The effect of the shape of six different periodic forces and second periodic forces on the onset of horseshoe chaos 
are studied both analytically and numerically in a Duffing oscillator. The external periodic forces considered are 
sine wave, square wave, symmetric saw-tooth wave, asymmetric saw-tooth wave, rectified sine wave, and modulus 
of sine wave. An analytical threshold condition for the onset of horseshoe chaos is obtained in the Duffing oscillator 
driven by various periodic forces using the Melnikov method. Melnikov threshold curve is drawn in a parameter 
space. For all the forces except modulus of sine wave, the onset of cross-well asymptotic chaos is observed just above 
the Melnikov threshold curve for onset of horseshoe chaos. For the modulus of sine wave long time transient motion 
followed by a periodic attractor is realized. The possibility of controlling of horseshoe and asymptotic chaos in the 
Duffing oscillator by an addition of second periodic force is then analyzed. Parametric regimes where suppression of 
horseshoe chaos occurs are predicted. Analytical prediction is demonstrated through direct numerical simulations. 
Starting from asymptotic chaos we show the recovery of periodic motion for a range of values of amplitude and 
phase of the second periodic force. Interestingly, suppression of chaos is found in the parametric regimes where the 
Melnikov function does not change sign. 
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1. Introduction 



The dynamics of nonlinear oscillators and circuits are very often studied with external force being of the 
form /sinwt or / cos tot. Other forms of forces such as square wave, rectified sine wave, rectangular wave, 
etc. can be generated and applied to dynamical systems. The study of the effect of such periodic forces will 
be helpful to choose a suitable external drive in creating and controlling nonlinear behaviours. In recent 
years there are reports on the effect of different forces on certain nonlinear phenomena [1 2 3 4 5 6J . Analysis 
of features of a particular dynamics with various periodic forces and a detailed comparative study of effects 
induced by them will be of great use. It is also important to explore the utility and applicability of analytical 
methods such as multiple-scale perturbation theory and Melnikov method to the systems driven by periodic 
forces other than /sinwi and / cos tot. 
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Motivated by the above, in the present paper we wish to study the occurrence of horseshoe chaos in 
Duffing oscillator driven by different periodic forces applying Melnikov analytical method. Recently, this 
method has been applied to study the Duffing and Duffing-van der Pol oscillators with different external 
perturbations [7}819|10|lljl2|13|l4]. The Melnikov method has also been applied to other nonlinear systems. 

In the present paper, we consider the perturbed Duffing double-well oscillator 

x = y , (la) 
y = U%x - (3x 3 + e [-ay + F(t)] , (lb) 

where a is the damping constant, Uq is the natural frequency, (3 is the stiffness constant which plays the role 
of nonlinear parameter, F(t) is an external periodic force, e is a small parameter, uiq > 0, (3 > 0. We study 
the occurrence of horseshoe chaos with different periodic forces such as sine wave, square wave, symmetric 
saw-tooth wave, asymmetric saw-tooth wave, rectified sine wave and modulus of sine wave. In section 2, we 
obtain the Melnikov threshold condition for the transverse intersection of homoclinic orbits for the system 
(p} separately for each of the above periodic forces. In section 3, we plot the Melnikov threshold curve in the 
(/ — ui) parameter space for all the forces where / and oj are the amplitude and frequency of the external 
periodic force. We verify the analytical prediction with the numerically calculated critical values of / at 
which transverse intersections of stable and unstable manifolds of the saddle occur. The Melnikov threshold 
value is also compared with the onset of asymptotic chaos wherever possible. Only for the modulus of sine 
wave long time transient motion instead of asymptotic chaos near the Melnikov threshold curve is found. We 
characterize the transient dynamics using the average transient time. In section 4, we consider the system 
((T|) driven by two sine and modulus of sine forces and analyze the effect of the amplitude and phase of the 
second force. Starting from horseshoe chaos and asymptotic chaos we show the possibility of suppressing 
them. The control regions can be identified by the Melnikov method. Finally, we end up with conclusion in 
section 5. 

2. Calculation of Melnikov function 

The unperturbed part of the system |T]) with e=0 has one saddle point (x*,y*) — (0,0), and two center 
type fixed points, (x*,y*) =(±a/wq//3, 0). The two homoclinic orbits connecting the saddle to itself are given 

by 



W ± {x h (t), y h {t)) = [±^2ujI/I3 sech 

=F w o sech [^f^o t) tanh ( \/ u o t ) ) ' ( 2 ) 

The Melnikov function M(to) measures the distance between the stable manifold (Ws) and unstable manifold 
(Wu) of a saddle. When the two orbits are always separated then the sign of M(to) always remains same. 
M(to) oscillates when the orbits Wu and Ws intersect transversely (horseshoe dynamics). A zero of M(to) 
corresponds to a tangential intersection. The occurrence of transverse intersections implies that the Poincare 
map of the system has the so-called horseshoe chaos [1115116] . 
For the Duffing eq.Q}, the Melnikov function is 

oo 

M[t Q )= J y h [-ay h + F(T + t )]dT. (3) 

— oo 

In the following we calculate the Melnikov function for the system ([T]) with different periodic forces. 
For the system ([T]) driven by the force F(t) — / s j n sin cot, the Melnikov integral §5§ is worked out as 

M L(to) = A± f sin B cos Lut Q , (4a) 
where 



2 



For the square wave F(t) = F S q(t + 2it/uj) = f S q sgn(sinwi) where sgn(y) is sign of y, 
series we obtain 

oo 

^sq(*o) = A ± / sq B n cos(2n - l)wt„, 
where 

4V2 w , / (2n- 1) 
B n = — -=— sech 1 ' 



For the symmetric saw-tooth wave of the form 

[ 4/sst^/r, < t < 7T/2LU 

F SS t(t) = F S st (t + V") = | -4/ sst i/T + 2/ sst , it/2lu < t < 3tt/2lu 

{ 4/sstVT - 4/sst, 3tt/2 W < t < 2tt/u> 

the Melnikov function is 



oo 



n=l 

where 



" n ttVP (2n-l) S6Ch 



(2n — 1)7tcj 



2v^ 



For the asymmetric saw-tooth wave 

^ast(*)= f, ast (* + 2tt/w) = 
we obtain 



2/ast*/T, < t < tt/w 

2/ast*/^ - 2/ast. < 1 < 



00 9/9 / 

M± t (t ) = A ± / ast ^ B n cosn^o, B« = -^=f sech " " 



\2^ 
For the rectified sine wave 



Free (t) = -Free (t + 2tt/w) = 

we find 

00 

MteSo) = A ± 2/ rcc ^ B n C„ sinnut , 

n=l 

where 

(2n-l)7r/oi 



/recsinwi, < £ < 7t/cj 
0, 7r/w < i < 27r/w 



VSlouq 1 

(2n-2)w/w 



B„ = _" — -. f'„ = y seel) \ \J ' J^t j laiil.) { \J J sin n^rdr 



For the modulus of sine wave -FmswW = ^msw(* + 2it/uj) = /msw sin (wt/2), we find 



oo 



^msw(*o) 



A ± /msw > _ -B„ sin nwio , 



(11a) 



n=l 



5, 




) 



(lib) 



3. Horseshoe chaos and strange attractor 

In this section, we compute the Melnikov threshold values for horseshoe chaos and compare it with 
numerical prediction. For the sine wave force the threshold condition for transverse intersections of stable 
manifolds (W^) and unstable manifolds (Wjf ) is 



For the other forces M(to) is a convergent series. In these cases the threshold values of / can be determined 
numerically. For the rectified sine wave force C„ 's in eq. (llObj) have to be calculated numerically. Throughout 
our study we fix the parameters in eq.JT]) as a — 0.5, f3 = 1, uiy = 1 and u = 1. /m can be calculated 
numerically as follows. For the forces other than sine and modulus of sine wave, we cannot write the 
sufficient condition for the existence of simple zeros of M(to). Therefore, we identify the occurrence of 
homoclinic bifurcations by numerically measuring the time tm elapsed between two successive zeros of 
M(to). tm is calculated for a range of amplitude of the external force. The value of / at which first time 
M (to) changes sign and thereby giving finite tm is the Melnikov threshold value for homoclinic bifurcation. 

Figure |T]) shows the plot of /m versus n, the number of terms in the summation, cq. (|5a[) . for the square 
wave force, Jm converges to a constant value with increase in n. For n > 10, the variation in Jm is negligible. 
Similar result is found for the other forces also. In our numerical calculation of Jm we fix n — 50. Figure 
([2]) shows the plot of the threshold curves for horseshoe chaos in the (/ — io) parameter plane. Below the 
threshold curve no transverse intersection of stable and unstable manifolds of the saddle occurs and above 
the threshold curve the transverse intersection of orbits of the saddle occurs. The threshold curves are 
nonintersecting. Smooth variation of /m is found for all the forces considered in our study. The variation of 
/m with uj is similar for all the periodic forces. Among the six forces, /m is maximum for the asymmetric 
saw-tooth wave and is minimum for the square wave. Thus, onset of horseshoe chaos can be either delayed 
or advanced by an appropriate choice of periodic force. 

We verify the analytical prediction by numerically computing the stable and unstable manifolds of the 
saddle. In fig. we plotted the orbits of the saddle for two values of / — one for / < Jm and another for 
/ > /m for each of the forces. For clarity only part of the manifolds are shown. In the left side subplots, 
for / = 0.2, the stable and unstable orbits are well separated. In the right side subplots, for / > /m, we 
can clearly notice transverse intersections of orbits at certain places. The critical values of Jm and /cwc (at 
which onset of cross- well chaos occurs) for various forces are given in Table 1. The numerical result agrees 
well with the theoretical prediction. 

For the system ([T]) driven by modulus of sine wave onset of cross- well chaos (asymptotic) is not observed 
near the Melnikov threshold value. However, transient motion followed by a long time periodic motion and 
hysteresis phenomenon are observed. Different paths are followed when / is increased from to 0.5 and 
decreased from 0.5 to 0. Onset of jumping motion from left-well to right-well is found to occur at /j. = 
0.440547. Long time transient motion is found for / values just below For a value of /, we calculated the 
transient time tt, the time taken to settle on the right- well attractor for a set of 900 initial conditions in 
the region x E [—0.9, —0.6], x E [—0.2, 0]. Then, average transient time r is calculated for a range of values 
of /. Figure (j4]) shows the variation of r as a function of (/ — /<•)• t is found to scale as 1.04166(/ — / c )~ ' 49 - 




(12) 
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Table 1 

Critical values of and /cwc (the values of / at which cross- well chaos occur) for the Duffing cq.([T} for various forces with 
a = 0.5, wjj = 1.0, = 1.0 and u = 1.0. 

No. Force f M /cwc 

1. Sine wave 0.3565 0.3833 

2. Square wave 0.2745 0.2910 

3. Symmetric saw-tooth wave 0.4629 0.4792 

4. Asymmetric saw-tooth wave 0.7155 0.7326 

5. Rectified sine wave 0.5026 0.5192 

6. Modulus of sine wave 0.3985 



4. Effect of second periodic forces 

In this section, we consider the problem of controlling of chaos by the addition of second weak periodic 
forces. We consider only the two forces, namely, sine wave and modulus of sine wave. The effect of other 
forces can also be studied. 

In eq.([T]) we choose F(t) — f s - ln sinwt + g s i n sm(ilt + cf>) where <fi is the phase shift. Recently, Zambrano et 
al [17] studied suppression of chaos by the phase difference between the two external periodic forces. Wang 
et al [TT] found the occurrence of strange non-chaotic attractor in the Duffing oscillator for a range of values 
of phase </> with the force / cos(wt) + gcos(Clt + (t>(t j). 

For the system (fTJ), the Melnikov function is 

M s1n^°) = A±B /sm coswio ± C g sin cos(fii + 4>), (13a) 



where A and B are given by eq. (|4b[) and 

C = Trfl sech ( -^L ) . (13b) 



The choice f2 = to and <fi = is trivial. For f2 ^ u> a Melnikov condition similar to eq. (fl~2|) cannot be written 
for arbitrary values of a; and O. So, we numerically compute the time tm and identify the parametric regime 
where 1/tm ~ 0. In the absence of the second periodic force the horseshoe chaos is found to occur for 
/=0.359 and uj = 1. We study the possibility of suppression of chaos by the addition of the second periodic 
force. Figure shows the plot of l/rj\/ against fl for g = 0.1 and cf> = 0. In this figure 1/tm is zero (that is, 
tm is infinity) only at f2 = 0.5, 1.5 and £1 > 2.0. This implies that horseshoe chaos does not occur for these 
three values of fL For other values of fi, both M + (£ ) and M~(t ) oscillate and hence 1/tm are nonzero. 
These results are verified numerically. 

Next, we study the effect of (j>. For u> = and ^ (which gives B — C), the condition for M(to) to 
change sign is given by 

ffs 2 in + 2/ sin 3sin cos + / s 2 in - (A 2 /B 2 ) > 0. (14) 

The equality sign corresponds to tangential intersections of the orbits of saddle. We denote g^ in and g~ in are 
the two roots of the eq. (fT4")) with equality sign. 

Figure ([6]) shows the plot of g s i n versus cf> for u> = O = 1, and / = 0.359. In the regions a, b and c enclosed 
by the continuous and dashed curves in fig.©, the sign of M(to) remains same and transverse intersection of 
orbits of the saddle does not occur. In the remaining region M(to) oscillates and horseshoe chaos occurs. The 
above result is also confirmed by direct numerical simulation. We note that suppression of horseshoe chaos 
can be achieved for a range of values of 0, f2 and g sm . A striking result is that suppression of asymptotic 
chaotic motion is also found in the regions a, b and c. In fig.® bifurcation diagram and maximal Lyapunov 
exponent (A) as a function of g s i n are reported. For 4> — 0.5 the maximal Lyapunov exponent A is found to 
be negative in the interval g s - ln £ [—0.6077, 0.0108] while the interval of g in which suppression of horseshoe 
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chaos predicted by the Melnikov method is [—0.6097,0.0103]. For <fi = it, suppression of asymptotic chaos 
is found in the interval g S i n G [—0.017,0.735]. The Melnikov method predicted interval of g for suppression 
of horseshoe chaos is [—0.015, 0.735]. Similar results are observed for various values of <fi m the regions a, b 
and c. 

For the Duffing oscillator driven by two modulus of sine wave 

oo oo 

A4 sw (t ) = A± /msw B n sm(ncut ) ± g m sw ^ Dn sm ( n ^ t o + <t>) , (15) 

n—l n—1 

where A and B n are given by the eqs. (|4bj) and (jTTJ) respectively and D n is same as B n except that u> 
in B n is now replaced by f2. We numerically compute the time tm and identify the parametric regions 
where 1/tm ~ 0. Figure ([5]) shows the plot of threshold values of gmsw versus <f> for u> = Q = 1, and 
/msw = 0.5. In the regions a, b and c enclosed by the continuous and dashed curves in fig.©, the sign 
of M(to) remains same and transverse intersection of orbits of the saddle does not occur. In the remaining 
region M(to) oscillates and horseshoe chaos occurs. The above result is also verified by numerical simulation. 
For example, fig. ([9]) shows the part of the stable and unstable orbits in the Poincare map for two values of 
<?msw chosen outside and inside the region b for u) = Q = l, 4> = tt and / = 0.5. The threshold values of 
<?msw are <? msw = 0.615 and g~ in = 0.393. When (?msw is varied from say —0.7 horseshoe chaos does not 
occur in the interval g m sw < 9 < ffmsw- I n n g-(EM transverse intersections are seen for g m sw = —0.1. For 
5msw = 0.5 (lying in the interval [g m sw; <?mswD t ne stable and unstable orbits are well separated (figOa). 
Suppression of horseshoe chaos can be achieved for a range of values of <f>, f2 and <?msw- The bifurcation 
diagram shows absence of cross-well chaos (asymptotic) near the Melnikov threshold value. However, long 
transient motion followed by a periodic motion and hysteresis phenomenon are observed. 

5. Conclusion 

In the present paper we have performed a numerical and analytical studies of homoclinic bifurcation in 
the Duffing oscillator driven by different periodic forces and addition of second periodic forces. Applying 
the Melnikov analytical method we obtained the threshold condition for onset of horseshoe chaos, that 
is, transverse intersections of stable and unstable branches of homoclinic orbits. For the modulus of sine 
wave instead of cross-well chaos a long time transient motion followed by a period-T solution is observed. 
Analytical prediction of horseshoe chaos is found to be in good agreement with numerical simulation for all 
the forces. Suppression of horseshoe chaos and asymptotic chaos are found for a range of amplitude of the 
second periodic force and phase shift between the two external forces. With the good agreement obtained 
between theoretical and numerical predictions we emphasize that the Melnikov analysis can be successfully 
used to predict the onset of chaos in the presence of weak periodic perturbations and also the effect of such 
perturbations on regular and chaotic dynamics. 

In the Duffing oscillator driven by each of the different periodic forces considered in our study transverse 
intersection of left-well homoclinic orbits and transverse intersection of right-well homoclinic orbits occur 
at same critical values of, say, the amplitude / of the force. This is because in the absence of driving force 
and damping the potential is symmetric and the two homoclinic orbits connecting saddle to itself satisfy 
the relation W~ = —W + . It is of interest to consider the situation where W~ ^ ±IF + . In this case, it is 
possible to have onset horseshoe chaos (as well as asymptotic chaos) in the two wells at different threshold 
values. This will be analyzed in future. 
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Fig. 1. /jvf versus n, the number of terms in the summation in eq, ]5aH . for a = 0.5, f3 = 1, Wg = 1 and w = 1 when the system 
JTT is driven by the square wave force. Variation in Jm converges to a constant value with increase in n. 
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Fig. 2. Melnikov threshold curves for horseshoe chaos in the (/ — ui) plane for the system JTJ driven by the forces (a) asymmetric 
saw-tooth wave, (b) rectified sine wave, (c) symmetric saw-tooth wave, (d) modulus of sine wave, (e) sine wave and (f) square 
wave. The values of the parameters in eq. JTJ are a = 0.5, /J=land^ = 1. 
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Fig. 3. Numerically computed stable and unstable manifolds of the saddle fixed point of the system JTJ. The system is driven 
by (a-b) sine wave, (c-d) square wave, (e-f) symmetric saw-tooth wave, (g-h) asymmetric saw-tooth wave, (i-j) rectified sine 
wave and (k-1) modulus of sine wave. The other parameters are a = 0.5, /3 = 1, u>^ = 1 and uj = 1. Left side subplots are for 
/ < /m while the right side subplots are for / > 
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Fig. 4. Variation of r as a function of (/ — / c ). r is found to scale as 1.04166 (/ — /e) — ' 
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Fig. 5. 1/tm versus f2 for the system JT} with -F(t) = f aln sinujt + g s i n sin(f2t + 0), g s in = 0.1, / s i n = 0.359, a = 0.5, j3 = 1, 
£Jq = 1, £J = 1 and = 0. Continuous curve corresponds to positive sign while dashed curve corresponds to negative sign in 
eq.|[T3t. 
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Fig. 6. Graph of the roots of g ain of ea. <|14|l with equality sign. Horseshoe dynamics does not occur in the regions a, b, 
and c enclosed by the curves. Continuous and dashed curves represent and g~ respectively. The other parameters are 
/ Bin = 0.359, a = 0.5, ft = 1, ui = fi = 1 and = 1. 
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Fig. 7. Bifurcation diagrams and the corresponding maximal Lyapunov exponent illustrating the effect of <j> and g on chaotic 
dynamics. The parameters are fixed as a = 0.5, /9 = 1, U)q = 1, U) = Q = 1, / s i n = 0.359. For <? s i n = the system exhibits 
chaotic motion. 
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Fig. 8. Graph of <? msw . Horseshoe dynamics docs not occur in the regions a, b, and c enclosed by the curves. The other 
parameters are /msw = 0.5, a = 0.5, f) = 1, ui = CI = 1.0 and uJq = 1. 
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Fig. 9. Numerically computed part of perturbed homoclinic orbits for (a) jmsw = —0.1 and (b) gmsw = 0.5. The other 
parameters are fixed at a = 0.5, (3 = 1, u)q = 1, (p = n, CI = w = 1 and /msw = 0.5. 
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